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The exclusive processes e^e^ VP, in the region of which the final state meson momentum 
is much larger than the hadronic scale Aqco, are studied in the framework of PQCD approach 
based on the kr factorization. Including the transverse momentum distribution in the light cone 
wave functions, our results are consistent with the experimental measurements. According to our 
results, many processes have large enough cross sections to be detected at ^ = 10.58 GeV. The 
s dependence of the cross section has been directly studied and our result indicates that the 1/s^ 
scaling is more favored than l/s'*. We also find that the gluonic contribution for the processes 
involving 77' ' is tiny. 

PACS numbers: 



I. INTRODUCTION 



In exclusive or inclusive processes with large momentum transfers, the production rates and many other 
phenomena, such as the dimensional rule, the hclicity structure, can be successfully explained by the pcrturbative 
QCD analysis . The essential ingredient is the factorization theorem which insures that a physical amplitude 
can be represented as a convolution of a hard scattering kernel and hadron distribution amplitudes. The former 
can be calculated using the perturbation theory while the latter, although non-pcrturbative in nature, are 
universal. The light cone distribution amplitudes which describe the longitudinal momentum distribution of 
partons in hadron, can be determined by the experiments of various channels. In e~^e^ — > 'y* VP at high 
energies (V denotes a light vector meson and P denotes a light pscudoscalar meson), the energy of the light 
meson is much larger than its mass and the hadronic scale Aqcd- One important feature of this process is that 
the meson moves nearly on the light cone. The energetic light meson is composed of two valence quarks which 
are both energetic and coUinear. The gluon which generates the quark pair is very hard and this leads to the 
application of pcrturbative QCD into this process. However, the application of the pcrturbative QCD approach 
to this simple process is complicated by the end point problem. If collincar factorization is applied, the hard 
kernel contains the inverse term of momentum fraction which makes the integration divergent at the end point. 
This divergence arises from the overlap of the soft and collincar momentum region ^. A modified pcrturbative 
QCD approach based on kx factorization, which keeps the intrinsic transverse momentum of partons in the 
meson, is proposed and successfully applied to many processes 0, H] ■ In this approach, the Sudakov effect is 
taken into account and the applicability of pcrturbative QCD can be extended down to a few GeV scale. It is 
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^ This overlap has also been attempted to subtract out in Ref. Q and new factorization theorems in rapidity space are subsequently 
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claimed that the perturbative calculations could be consistent at scale about Q ~ 20Aqc'd in this framework. 
This approach is also called PQCD approach for simplicity. 

The exclusive two-meson productions in e^e^ annihilation provide an opportunity to investigate the behaviors 
of various meson form factors. The dependence of the form factor on the energy scale can shed light on the 
internal strong interaction information. It can also give information on the wave function of the hadrons in 
terms of its partonic constituents. In the standard model, the exclusive production of hadron pairs at e^e~ 
colliders can proceed through a virtual photon or a boson. At energies well below the mass of Z", the 
production proceeds predominantly via the annihilation of e+e^ into a virtual photon. Due to the invariance 
of charge conjugation in electromagnetic and strong interactions, the final state should have the same charge 
conjugation quantum number with a photon, i.e., these processes can only produce final states with charge 
conjugation quantum number C = — 1. e+e^ VP can proceed via the following form factor: 

(F(6,pi)P(p2)|jriO) = Fvp{s)e^.o.pe'^p'ipl (1) 

where Pi{p2) is the momentum of the vector (pseudoscalar) meson and e is the polarization vector of the vector 
meson. Here j^™ is defined as j^™ = 97/^9- Eq. ([T]) indicates that the vector meson is transversely polarized. On 
the experimental side, the productions of VP have been extensively studied: BES and CLEO-c have reported 
the continuum productions 3, Q, Recently, BaBar Collaboration observed the exclusive reaction e+e" — > (firj 
at ^/s = 10.58GeV and measured the cross section Since fc^ factorization can give a reliable prediction in 
other similar processes, in this paper, we will perform a study on e+e^ —t^*^ VP in this framework and 
make a comparison with the data. 

Another interesting reason for investigating the e+e~ annihilation is the similarity with the annihilation 
corrections in B decays, shown in Fig. [TJ In charmless two body B decays, the annihilation diagrams are 
power suppressed relative to the emission contribution. But it is found that in B —>■ ttK, tttt decays, the 
factorizable annihilation diagrams could be important due to the chiral enhancement factor for operator Oq . 
This enhancement of factorizable diagrams can provide large strong phases and give large CP asymmetries, 
which indicates the annihilation diagrams are of great importance. Comparing the two diagrams in Fig. [H we 
can see that the e~^e~ annihilations have similar topologies with the factorizable annihilation diagrams in B 
decays. They may provide an ideal laboratory to isolate the power correction effect and to find out whether 
contributions of annihilations from end-point are important or not for meson productions [lol | . 






(a) 



FIG. 1: The annihilation diagrams in B decays and e+e~ annihilation. In the left diagram, B meson is annihilated 
through the four-quark operator. In the right diagram, electron and positron annihilate into a virtual photon. These 
two diagrams have similar topologies. 
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TABLE I: Input values of the decay constants of the pseudoscalar and vector mesons (in MeV) 






fK 


fp 


Jp 






/a- 


Ik* ftp 


/J 


131 


160 


209 


165 ±g 


195 


145 ± 10 


217 


185 ± 10 231 


200 ± 10 



The remainder of this paper is organized as foUows. In Section [TTl we present the expression for the cross 
sections for e+e^ VP in the kx factorization: the first part of this section is devoted to the discussion on the 
decay constants and the distribution amphtudes of mesons, and the second part contributes to a brief introduc- 
tion to PQCD approach and factorization formulae for form factors. The numerical results and discussions are 
presented in Section [IIII The last section is our summary. 



II. CALCULATION IN kr FACTORIZATION 



A. Decay constants and Wave functions 



The decay constants for a pseudoscalar meson and a vector meson are defined by: 

1927^.759110) = -*/pP^, (l/(P,e)|g27M9i|0) = /ymye^, {V{P,e)\q2<7^M^) = -^fvi^f^P^ ~ e-^^). (2) 

The pseudoscalar decay constants taken from the Particle Data Group [ll| are shown in Table HI The charged 
vector meson longitudinal decay constants are extracted from the data on —>■ {p^ , K*^)vt, while the neutral 
vector meson longitudinal decay constants are determined from the data on the electromagnetic annihilation 
processes — > e+e" 11 1. The transverse decay constants are taken from the QCD sum rules [3, which 
are also collected in Table HI 

The light-cone distribution amplitudes are defined by the matrix elements of the non-local operators at the 
light-like separations with = 0, and sandwiched between the vacuum and the meson state. The two- 
particle light-cone distribution amplitudes of an outgoing pseudoscalar meson P, up to twist-3 accuracy, are 
defined by Q: 



75 



^(x) -I- mo75(/>"^(a;) - moa^''j5P^z^ 



Q/3 



a0 



(3) 



where n, v are two light-cone vectors. The pseudoscalar meson is moving on the direction of n, with v the 



opposite direction, mo 



Mi 



is the chiral enhancement parameter, x is the momentum fraction carried by 



the positive quark q2. We have performed the integration by parts for the third term and <j>^{x) — ^-^(jf [x). 
The explicit form of distribution amplitudes for pseudoscalar mesons have been studied in QCD sum rule 
approach and other methods Q- I'^ principle, they are factorization scale dependent. Here we use the 
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following form for leading twist distribution amplitudes 



cj>^{x) = ^x(l-x)[l+ajcf^(t)], (4) 

cpiix) = ^x{l^:,)[l + a^cl'\t)+4cl'\t% (5) 

where t = 2x — 1 and Gegenbauer polynomials defined as: 

Cl'\i)^?,t,Cl'\t) = l^{bt^-l). (6) 

The Gegenbauer moments at /i = 1 GcV arc determined as 

flj'^ = 0.25 ±0.15, af = 0.06 ±0.03. (7) 

Since the momentum transfer at ^/s = 10.58 GeV is large enough, the use of asymptotic forms for twist-3 

distribution amplitudes is acceptable. Besides, we also use these forms at a/s = 3.67 GeV for simplicity. The 
asymptotic forms of twist-3 distribution amplitudes are given as: 



As for the mixing of 77 and 77', we use the quark flavor basis proposed by Fcldmann and KroU 17|, i.e. these 
two mesons are made of nn = {uu + dd)/\/2 and ss: 





= U{6) , (9) 

with the mixing matrix, 

/ cos — sin \ 

(10) 




where the mixing angle 9 = 39.3° ± 1.0°. In principle, this mixing mechanism is equivalent to the singlet and 
octet formalism, which is shown in [l^ . But the advantage is transparent, since only two decay constants are 
needed: 

{OWlBn\Vn{P)) = ^fnP^, 
{0\sri5s\Vs{P)) = ifsP^ . (11) 

We assume that the wave function of nn and ss is the same as the pion's wave function, except for the 
different decay constants and the chiral scale parameters: 

/„ = (1.07 ± 0.02)A, = (1.34 ± 0.06)/,. (12) 

The chiral enhancement factors are chosen as 

= [inl cos^ 9 + rn^, sin^ 9 - {m^, ~ ) cos 9 sin 6*] , (13) 

zniji jfi 

mff = —\ml, cos^ 9 + ml sin^ 9 ^-^{ral, - ml) cos6lsin6'l, (14) 



with m„ 5.6 McV and = 137 McV at ^ = 1 GeV [3- 

In this work, we also investigate thcgluonic contribution for iso-singlct pscudoscalar meson rj and rj' . This 
contribution has been attempted in |18j with neghgible effect in B rj form factor and a few percents to 
B ^ rj' . The leading- twist gluonic distribution amplitudes of the 77,1 and 77s mesons are defined as [l9|: 



MP)\Al{z)Al^{Q)\0) 



V3 AViN^-l 
fs Cf 



n- ■ P 



x{\ — x) 



V34V3iV2-l 



n_ • P 



x{l — x) ' 



where A 



''^^{z)Al^{w) = [Al{z)Al{w) - Al{z)A^^{w)]/2 and the function [20| 



= x\l - xfB-^'^Cl'\2x - 1) , Cl'\t) = 5t 



'5/2, 



(15) 



(16) 



The gluon labeled by the subscript ^ carries the momentum fractions x based on the above definition. The two 
Gegenbauer coefficients B2 and i3| could not be the same in principle. However, it is acceptable to assume 
B2 = i3| = -82, since there are large uncertainties in their values. Here the range of B2 has been extracted as 
B2 =4.6 ±2.5 Q. 

Following the similar procedures as for the pscudoscalar mesons, we can derive the vector meson distribution 
amplitudes for the transverse polarization up to twist-3 2l|: 

{V[P,eT)\q2p{z)qio.{m) = j\xe"''-^ [My Mv{x)+ /t ^vi^) 

+ Mvie^upal5l''^Tn''v''4>l{x)]^^ , 

where we have adopted the convention e^i^a _ ^ for the Levi-Civita tensor . 

The twist-2 distribution amplitudes for transversely polarized vector can be expanded as: 

fT 

V6- 
3/J 



(17) 



%x{l-x) \i + 4,cl'\t) 



c^l,{x) 



^x(l-x-) [l + a^^Cl'\t) 
^.x(l - ^) [1 + ai,,,Cl'\t) + 4K.Cl'\t) 
^^L(l-x-) \l + a^^Cl'\t) 



V6' 



The Gegenbauer moments have been studied extensively in the literature 
updated form Q, Q : 



21 



(18) 
(19) 
(20) 
(21) 

, here we adopt the very recent 



ai^, = 0.04 ± 0.03, a 



2p ~ "■2uj 



aL = 0.15 ± 0.07, a^i^. = 0.11 ± 0.09, a^^ = OMto fj- 



(22) 



As for the twist-3 distribution amplitudes (jiy and 0y, there is no recent update associate with those updates 
for twist-2 distribution amplitudes [l^. [isll. we also use the asymptotic form: 



1 

8V6 



[\ + {2x-l)% 0^(x) 



3/y 
4^/6 



(1 - 2x). 



(23) 
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The above discussions concentrated on the longitudinal momentum distribution and we intend to include 
the transverse momentum distribution functions of the pseudoscalar and vector mesons. But at present, the 
intrinsic transverse momentum dependence of wave function is still unknown from the first principle of QCD. 
As an illustration, we use a simple model in which the dependence of the wave function on the longitudinal and 
transverse momentum can be factorized into two parts 



231: 



ip{x,'kT) = (/>(x) X S(kT), 



(24) 



where 4'{x) is the longitudinal momentum distribution amplitude which has been discussed above and S(kT) 
describes the transverse momentum distribution. S(kT) satisfy the normalization conditions: 



In the following, we will use a Gaussian distribution: 



E(kT) - ^exp(-/3^k^), 



(25) 



(26) 



where the parameter (3 characterizes the shape of the transverse momentum distribution. The numerical value 
for /3 can be fixed by the condition that the root mean square transverse momentum (k^)^/^ should be at the 
order of Aqcd- Their relation can be derived from: 



(k|) 



/„'rfx/d^kTk||V;(x,kT)P _ J_ 



(27) 



If we choose the root mean square transverse momentum (k^)^/^ = 0.35 GeV, then = 4GeV~^. In PQCD 
approach, the integration will be transformed to the b space (coordinate space) and it is convenient to use the 
Fourier transformation of S(k7-): 

f,2 



S(b)= J d2kTe-'^-^E(kT)=exp(-A_). 
It can be observed that in the limit /? — > oo, S(b) can be simply replaced by 1. 



(28) 



B. Form factor and cross section in fcy factorization 



In the center of mass frame, we define qi, q2, pi and p2 to be the four- momenta of e^, in initial states, 
vector (V) and pseudoscalar meson (P) in final states, and define fci(2) and Xi(2) to be the momenta and 
momentum fractions of the positive quarks inside V and P respectively. The center mass energy of this process 
is denoted as Q = ^/s. Using the definition of the form factor in Eq.(l), we can obtain the cross section as 



VP) 



6 



IF, 



VP\ 



2<i>3/2(s), 



with 



$(,s) = 



1 - 



{mv + mpY 



1 - 



(mv — m,pY 



(29) 



(30) 
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There arc four different types of diagrams contributing to the production of vector and pseudoscalar meson in 
e~^e~ annihilations, to the leading order of the strong and electromagnetic coupling constants. The first type of 
diagrams contributing to this process are displayed in Fig. [21 These diagrams give the dominant contribution. 
The diagrams in Fig. [3] contribute to the processes involving 77 and 77', while the diagrams in Fig. 3] only 
contribute to the processes involving a neutral vector meson, such as p". uj and (f). Although these diagrams 
are suppressed by Uem, they can be enhanced by s/Ag^^. This mechanism is similar with the enhancement 
in penguin-dominated B decays 



24| and the so-called fragmentation mechanism in e+e VV processes 



25 



It is also interesting to explore this effect in e+e 



VP. For e+e~ 



K*K and e+e' 



the two 



photon non-fragmentation diagrams can give their contributions as in Fig. [5l But these diagrams suffering the 
suppression from electromagnetic coupling constant a^m which can be neglected safely. 

V 




(a) 




(b) 




(c) 




(d) 



V 



FIG. 2: Dominant contribution of e e 



VP 



We begin with a brief review of the PQCD approach. The basic idea of PQCD approach is that it takes 
into account the transverse momentum of valence quarks which results in the Sudakov factor. The form factor, 
taking the first diagram in Fig. [5] as an example, can be expressed as the convolution of the wave functions tpvj 
ipp and the hard scattering kernel Th by both the longitudinal and the transverse momenta: 



Fa{VP)= j dxidx2 J d^kTld^kT2^pv{xl,'kTl,Pl,^J.)TH{xl,X2,Q,kTl,'kT2,^^)i^p{x2,'kT2,P2,^J')■ 



ThIO\lg]l the Fourier transformation, the above equation can be expressed as: 



Pais) 



dxidx2 



;Vv{xi,hi,pi,fl)THixi,X2, (3,bi,b2,/i)Pp(x2,b2,_P2,M)- 



(31) 



(32) 



Here Vi{xj,hj,pj, 11) are the Fourier transformation of il}i{xj,'k.TjTPj, jj), where the subscript i denotes V or P, 
and j indicates 1 or 2. 

In the above expression, the double logarithms, ari sing from the overlap of the soft and coUincar divergence, 
have been resummed to result in the Sudakov factor 12611 



'Pi{xj,hj,pj,fi) = exp [-s{xj,bj, Q) - s(l - Xj,b.j, Q)] V^{xJ,hJ, fi) 



(33) 
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The exponent s(^, bj,Q), ^ — Xj or 1 — Xj, is expressed as 



In 



\V2pJ 



'i/b, P 

where the anomalous dimensions A and B to one loop are given by 

rv 2 (T /p27E-i 

A ^ Cf — , B = -^\n' 

TT 3 TT 

with Cp = 2N ^^'^ being the Eulcr constant. The one-loop running coupling constant, 



/3oln(/iVA2 



QCDl 



with the coefficients 



/3o 



33- 2n/ 



(34) 



(35) 



(36) 



(37) 



where Uf is the number of the active quark number. We require the relation of the involved scales S_Q/^/2 > 
1/bj > A as indicated by the bounds of the variable p in Eq. ((34|) . The QCD dynamics below 1/bj scale is 
regarded as being nonperturbative which can be absorbed into the initial condition Vi {xj , hj , ^) . 

The form factor, as a physical observable, is independent of rcnormalization scale /i, but the functions V and 
Th still contain single logarithms from ultraviolet divergences, which can be summed using the rcnormalization 
group equation method. This rcnormalization group analysis applied to Th gives 



THixj,hj,Q,fi) = exp 



^7g(as(M)) 



X THixj,hj,Q,t) , 



(38) 



where = — as/Tr is the quark anomalous dimension in axial gauge and t is the largest mass scale involved in 
the hard scattering, 



t = ma.x{y/x^Q, 1/bi, 62) 



(39) 



The scale ^Jx^Q is associated with the longitudinal momentum of the quark propagator and with the 
transverse momentum. The large-6j behavior of V is summarized as 



s{xj,bj,Q) - s{l - Xj,bj,Q) - 2 / — 7g(as(Ai)) 

1/6, M 



Viixj,hj,pj,fi) = exp 

where Vi{xj, hj, 1/bj) is the wave function discussed above: 

n(xj,b„ 1/5,) ^ Hx,,l/b,) X S(b,). 



The threshold resummation 



23, 



28j can also play an important role in 



x7^,(x,b,,l/6,), (40) 



(41) 



VP processes. The lowest- 



order diagrams Fig. [2lja) and (d) give an amplitude proportional to 1/(0:2(1 — xi)) and 1/(2:2X1) respectively. 
In the threshold region with X2 — > (to be precise, X2 ^ 0(Aqp^/s)), additional coUinear divergences are 
associated with the internal quark. The QCD loop correction to the electromagnetic vertex can produce the 
double logarithm as In^ X2 and resummation of this type of double logarithms lead to the Sudakov factor 5*4(2:2). 



Similarly, resummation of In^ xi due to loop corrections in the other diagrams lead to the Sudakov factor 
St{xi). The Sudakov factor from threshold resummation is universal, independent of flavors of internal quarks, 
twists, and the specific processes. To simplify the analysis, the following parametrization has been used 



with the parameter c ~ 0.3. This parametrization. symmetric under the interchange of x and 1 — x, is convenient 
for evaluation of the amplitudes. It is obvious that the threshold resummation modifies the end-point behavior 
of the meson distribution amplitudes, rendering them vanish faster at a: ^ 0. 

Combing all the above ingredients, we obtain the factorization formula for the contribution from Fig.[2]^a): 

/•l poo 

Fa{VP) = l&irCFQri / dxidx2 / 5id6i&2d&2 
Jo JO 

y«pj,{x2M)Wv{xi,bi) - 4>l{xiM)]E{ta)h{l - .ti, X2, &i, 62), (43) 

where h and E are defined by 



h{x^,X2MM) = C-^?St{x2)[e{h^-b2)H^^\^QbM^/^Qb2) 

+9{b2 - bi)H^o'\V^Qb2)MV^Qbi)\H^o'\V^Qbi), (44) 

E{ta) = as{ta)exp[~Si{ta)~ S2{ta)], (45) 

where Jq and H^^^^ are the Bessel functions, respectively, ta = max(y/2?2Q, 1/&2) and ri = My/Q- 
Similarly, for the other diagrams, the amplitudes are: 

poo 

Fb{VP) ^ -WttCfQ dxidx2 bxdbxb2db2E{tb)h{x2,l - xiMM) 
Jo Jo 

x{ri(xi-l)[(/)5'(xi,5i) + 05:(xi,6i)]0^(a;2,62) + 2r20f(xi,6i)(/)^(x2,62)}, (46) 

Fc{VP) = -IGnCpQ [ dxidx2 [ bidbib2db2E{U)h{\ - X2,xiMM) 
Jo Jo 

X [ria;i(0?(xi, 61) - 0i"(xi, 6i))0^(x2, 62) + 2r20f (xi, 6i)</)^(a;2, 62)], (47) 

n 1 nOC 

Fd{VP) = -WTTCpQn / dxidx2 / b,db,b2db2E{td)h{xi,l - X2,bi,b2) 

Jo Jo 

^{<t>l{x,) + cj>l{x^)^^{x2M), (48) 

with r2 — m^lQ. The factorization scales ti are chosen as 

U = max(\/l - xiQ, 1/62), tc = max(^i?7(3, 1/62), td = max(x/l - X2Q, I/&1, 1/&2). (49) 

If the final state meson is not K* or K, the distribution amplitudes are completely symmetric or antisymmetric 
under the interchange of Xj and I ~ Xj. Then one can easily obtain: 

Fa{VP) = FdiVP), Ft,{VP)=F,{VP). (50) 

For the flavor-singlet pseudoscalar meson ry and r;', there are additional contributions from the two-gluon 
diagrams as displayed in Fig. [31 even if they may be suppressed by the gluonic distribution amplitudes. However, 



10 




(a) 




(b) 




(c) 



FIG. 3: Gluonic contributions 



it is still worthwhile to investigate the numerical contribution in order to make our calculations as complete 
as possible. The computations of these diagrams are much similar with that showed in Fig. [2l The explicit 
calculations show that the Fig. [3] (a) docs not contribute to the transition amplitude, due to the antisymmetry 
of the two gluons. The amplitudes of the other two diagrams are given as 



X [(.T2 + 1)</)^(X1, 5i) ~ {X2 - 1)</)^(X1, 



— I dxidx2 / bidbib2db2E{te)h{x2,Xi,bi,b2) 
Ij Jo Jo 

X2[l - X2) ' 

— / dxidx2 / bidbib2db2E{tf)h{l - X2,l - xi,bi,b2) 
1) Jo Jo 

4'^{x2,b2) 



X (xi - l)[(a:2 - 2)(f>'"y{xi,bi) - a;2</)^(xi, 61)] 



3:2(1 - X2) ' 



for e+e Vrjg process with 



te = max(yxrQ, 1/61, 1/&2), tf = max(\/l - xiQ, l/bi, 1/62)- 



(51) 



(52) 



(53) 



It should be pointed out that the factor "2" from the exchange of two identical gluons in the final states has 
been added in the above equations. The amplitude for e+e~ — ^ Vrjn can be easily obtained by replacing the 
corresponding decay constant with an additional factor \/2 from Eq. (|51l52p . 

Furthermore, there are also contributions from the transition of photon radiated from one valence quark in 
pseudoscalar meson into vector meson directly, which have been presented in Fig. [4l Although these diagrams 
maybe suppressed by coupling constant of electromagnetic interactions, they are also enhanced by the almost 
on-shell photon propagator compared with the first type diagrams, especially for the processes with a very large 
center mass energy. These diagrams can also be calculated according to kx factorization, however, we will 
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{a) (b) 

FIG. 4; Enhanced diagrams for a neutral vector meson production 

simply adopt coUinear factorization due to disappearance of infrared divergence for these two diagrams. These 
two amplitudes are equal after integrating the momentum fractions carried by the valence quark of the meson. 
Hence, we obtain the amplitudes corresponding to them as follows: 




FIG. 5: Two photon non-fragmentation diagram. This contribution is suppressed by a^m and can be neglected. 
The form factors for the explicit channels can be easily obtained from the combinations of the eight amplitudes 



Fa-h- To be more specific, we can write them as 

i^p+x- = i[F,(p7r) + Fb(H], (55) 

Fp«.o ^ ^[Fa{pTi) + Ft,{p^)] + ^[Fg{pn)+Fu{pn)l (56) 

F^^o = [Fa{ujn)+F^{un)] + ^[Fg{pTT) + Fh{pTi)l (57) 

= -t|[^9(</'^)+^'^(<^'^)]' (58) 



Fk'^k^ = '^-[Fa{K*K) + F,{K*K)]-^-[F,{K*K) + Fd{K*K)l (59) 
F^.c^o - ^^[Fa{K*K) + F,{K*K)]-^[F,{K*K)+Fa{K*K)]. (60) 
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The form factor of e+e pjf- ^ can be written as the combination of its nn and ss component: 

Fy^ = COS0Fvrf„ - smOFy,,^, (61) 

Fvjj' = sinOFvrj^ + cos0Fvrj, , (62) 

where V — ,uj,(j) and 

FpO,,^ = [FaipVn) + FbipVn)] + ^[Fe(p?M) + Ff{p7]n)] + ^[FgipVu) + Fh{pTJ„)], (63) 

1 /2 
F,o,, = —[F,ipr,s) + Ffip^,)]-^[Fg{pT^,) + Fh{pVs)], (64) 

1 \/~^ 5 

^^r,„ = -[Fa(wr?„) + F,,(tJ7?„)] + -^[Fe(t^?y„) + i^/(w77„)] + — [Fg(t^7y„) + F?,(W7;„)], (65) 

/2 /2 
= ^[Fe(a.r,,) + F/(c^7,,)]-^[F<,(a.r,,,) + F^(c^r/s)], (66) 

F^r,^ = -i[Fe(077„) + F/(<^r;„)]-^[F<,(077„) + i^^(^77„)], (67) 

F^,u = + ^^6(<^^s)] - \\FM^s) + i^/(0??.)] - ^{F.A'i^^s) + F/.(</'?7.s)]. (68) 

III. NUMERICAL RESULTS AND DISCUSSIONS 
A. Cross section 



Making use of the distribution ampUtudes and the inputs listed before, one can easily obtain the cross sections 
for the process e+e^ PV. Here we would like to present the results of cross sections at =3.67 GeV and 
10.58 GeV in Table together with the data measured by CLEO-c and BaBar collaboration. The different 
scenarios SI, S2, S3 denoting different transverse momentum distribution functions, which will be discussed 
in the next subsection. As the longitudinal decay constants of the vector mesons and the pseudoscalar meson 
decay constants arc precisely determined, the uncertainties from these inputs arc neglected. Therefore the 
uncertainties shown in Table |TT] are from the transverse decay constants of the vector meson shown in table HI 

From the Eqs. ([55|) and ((56|) . we can see that, if neglecting the fragmentation contribution Fg^h, the cross 
sections for production of p+7r^ and p°7r'^ in e+e^ annihilation should be the same. At y/s ~ 3.67 GeV, the 
fragmentation can not give large contribution as the on-shellness enhancement is not strong. Thus theoretical 
calculation predicts that the ratio Ri = '^^^g+f.-'^piy^O) should be around 1 . From table |TT1 one can see that 
this prediction is consistent with the CLEO-c results. At higher energies, the enhancement effect becomes 
more important. This effect can weaken e^e^ — > p°7r° by about ten percents at ^/s = 10.58 GeV, relative to 
e^e^ —>■ p'^TT^. If the center mass energy is large enough, the contribution from diagrams in Fig. [4] will be 
dominant over the other contributions. 

The process e+e^ K*K has previously been calculated in PQCD {kx factorization) and has been shown 
to give correct order of magnitude for the form factors [2^. But they assume SU(3) symmetry using asymptotic 
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TABLE 11: Results of e^e~ VP cross sections at ^ — 3.67 GeV and y's = 10.58 GeV using three different transverse 
momentum distribution functions, denoted as SI, S2 and S3 respectively. Tlie experimental results from CLEO-c and 
BaBar collaborations are also included. 





= 3.67 GeV 


^fs = 10.58 GeV 


Channel 


(TQ1 fob) 


(Tc;'9 fob) 


f7<?QfDbl 
^ O O V_jr^ '-'_/ 


cjp-PT, fob) 


crsi(fb) 


(joo f fb) 


CT'^'? ffbl 

^ o o \ / 


I fb 1 


r " 


■J-"— 0.2 


-'-•■^-0.2 




^•"-1.2-0.5 


"• ' ^-0.04 


55+0'"? 


68+'!-0'; 

"■"°-0.04 




p " 


3 8+?, ? 

0.2 






q 1+1.0 + 0.4 

■'•^-0.8-0.4 


64+?,'nl 


50+?,;?:^ 


62+2 Si 






OQ + 2.2 
■^'-'•■^ — 2.2 


13 8+v; 


21 2+M 

— 1.6 


-1 r 9 + 2.8+1.5 


c 9+0.4 
■-^ — 0.3 


4 1+^1 

^•^ — 0.3 


5 0+°1 

^^ — 0.3 






1.2 X 10"* 


1.2 X 10~* 


1.2 X 10~* 


< 2.2 


2.1 X 10-3 


2.1 X 10-3 


2.1 X 10-3 




K*+K- 


5.6ir4 


2 q+°-i 

^•''-O.S 


4.31^3 


-1 n+i-i+o-5 

-'-•"-0-7-0.5 


-■ 9+0.02 
-'-•^-0.02 


0.83lH^ 








34.8t|j 


17.3ti;5 


26.4ti;| 


oq c+4. 6+3.1 
^■^•"-3. 9-3.1 


7 1+0.4 








P°V 


16.6t?;;; 


« 1+0-5 


12.5^0 7 


1"-U_i g_i 


q q+0.2 
'J"5-0.2 


2.41°;? 


q 1+0.2 
'^•-'--0.2 








*"J-0.3 


"•"-0.4 


9 -1+4.7+0.2 
^•^-1.6-0.2 


2 l+^'i 


1 5+°-i 


2 0+0-1 

^•"-0.1 




bjri 


1 5+0-1 


'"-0.07 


J^-J^-0.1 


9 Q+1.8+0.5 
^"^-1.0-0.5 


"•'5^-0.02 


r, 99+0.02 
"•^^-0.01 


29+0-02 

"•^■^-0.01 




ujrj' 


0.79iro^6 


"■•jy-0.03 


o.6ol°:°l 


< 17.1 


"•^"-0.02 


14+001 
"•J^^-0.01 


,-, -1 O+0.012 




(pV 


i9.ili;i 




14.6ir9 


9 1+1.9+0.2 
^- J^-1.2-0.2 


4 O+0.2 


q q+0.2 
0.0_o 2 


4 1+0-1 

^•J^-0.2 


9 Q+0. 5+0.1 
^•^-0.5-0.1 


cj>Ti' 


22.6tij 


11.51°;? 


17.4ll;J 


< 12.6 


r O + 0.3 
O-°-0.3 




c 4+0-4 





■wave functions. In order to sho-w the S'[/(3) symmetry breaking effect in e+e K* K , wc define the ratio: 

1 _L Fc+Fd 2 



(69) 

If "we assume that SU (3) symmetry "works -well, then the light cone distribution amplitude of K and K* is 
completely symmetric under the exchange of the momentum fractions of quark and anti-quark. We "will have 
Fa + Fb = Fc + Fd, then i?2 = 4 can be derived directly from Eq. ([59]) and One of the SU(3) symmetry 

breaking effects is that the s quark is heavier than n{— u,d) quark and carries more momentum in the final 
state light X^*-' meson. The gluon "which generates ss is harder than the nn generator, then the former coupling 
constant is smaller due to the more off-shell gluon. Consequently this leads to a smaller contribution to the 
form factor \Fa + Ft,\ than \Fc + Fd\- Therefore R2 is larger than 4. Using the cross sections listed in Table HIl 
we obtain our result for i?2: 

R2 = 6.0, (70) 

"where only the central value is given. The CLEO-c results indicate that there is a large deviation from the 
SU{3) limit J 

i?2 = 23.5^26.1 ± 12-2. (71) 

The central value of the experimental results for R2 seems too large, but as the uncertainties are also large, our 
result could be consistent "with results from CLEO-c collaboration. 

For the processes involving ry' such as the process e+e^ (f)!], wc find that the gluonic contribution is 
around one percent to the total cross section at ^/s = 10.58 GeV. This conclusion is consistent -with the study 
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on the B ^ T]^^ form factor [l^. The cross sections of e+e^ p{uj)ri'^'^ and lott^ at 3.67 GeV calculated 

in kT factorization arc consistent with the experimental values. The result for e^e^ (prj at ^/s ~ 10.58 GeV 
is also consistent with experimental data. This indicates that factorization is an effective method to deal 
with the infrared divergences in exclusive processes. 

The measurement of cross sections at different center mass energy -y/s can shed light on the s dependence. 
This dependence is expected as 1/ s3 or l/s^Q 

Isil l . Our results displayed in Table at two different scales 
y^s = 10.58 GeV and 3.67 GeV seem to favor the \/s^ scaling. It should be noticed that we neglected the 
Q ~ y/s dependence of the light cone wave function and ncxt-to-leading order contributions in our calculation. 
Therefore the s dependence study of the cross sections is not a complete one. 

As the quark or the gluon could be on-shell, we expect the amplitudes receive an imaginary part which is 



similar with the exclusive B decays PjIs^. The imaginary part in e^e~ (pij at y/s = 10.58 GeV is about twice 
as large as the real part in magnitude and a large strong phase is consequently generated. The contributions 
from Fig. [3] arc small; Contributions from Fig.[4]are small and real; The four diagrams in Fig.[2]give comparable 
contributions which are the main origin of the imaginary part. Unlike the B decays, strong phase here does not 
make any physical meaning, since there is no electroweak phase for interference. 

From Table ITU we can see that at y/s = 10.58 GeV, the cross sections for many processes, especially e+e" 
K^*K'^ and e+e" ^77', arc large enough to be detected. We suggest the experimentalists to measure these 
channels. 

In the above, we only concentrate on the exclusive production of a vector and a pseudoscalar. Applications to 
PP and VV productions are straightforward. The diagrams in Fig. [2] will give dominant contributions, where 
the final state must have negative charge conjugation quantum number C = — 1. Then only three channels for 
PP is allowed through one photon annihilation: e+e^ tt+tt^, e+e^ K^K^ and e+e^ K'^K'^ . If ?7-spin 
is well respected, d and s quarks are symmetric in K''^ and the cross section of e'^e^ K^K^ is zero. The 
non-zero result for e+e^ K^K^ can reflect the size of J7-spin symmetry breaking. For production of VV , 
the analysis is similar. Two flavor singlet vector mesons can not be produced through one photon annihilation 
diagrams either, but these productions could receive large additional contributions [25|. 



B. Theoretical Uncertainties 

One of the major uncertainties in our computations comes from the distribution amplitudes for the pseu- 
doscalar and vector mesons. The dependence on the longitudinal distribution amplitudes has been studied 



intensively in the exclusive B decays [33|. They will give 10-20% uncertainties here too. In the following, we 
will focus on the transverse momentum distribution. In PQCD approach, the intrinsic transverse momentum is 
taken into account. The resummation of large double logarithms results in the Sudakov factor which suppresses 
the large h region's contribution. As we can see from Eq. (|28p . the transverse momentum distribution hmction 
also suppress the contribution from the large b region. For momentum transfer of a few GeV, the transverse 



momentum distribution function damps more than the Sudakov factor 



23|. This suppression makes PQCD 



approach more self-consistent. So we expect that there is an obvious suppression for the production rate of 



15 



e^e^ — > VP at ^/s = 3.67 GcV if the transverse momentum distribution amplitude is taken into aceount. 
However, at present, it is still lack of first-principle study on the intrinsic transverse momentum distribution. 
The simple form is chosen as the Gaussian form discussed in Eq. (^5)1 or the following one, 



E(x, b) = exp 



x{l — x)b 



(72) 



4a2 

where a is the transverse size parameter as f3. As a simple test, we can choose a ~ I which is consistent with 



the value used in 



33|. Comparing with the form in Eq. ^E^, we can see that at x = 1/2, the two different 



forms coincides. But for small or large x, the second form can not give the same strong suppression as the first 
form (Eq. ([28]) ). We can expect the suppression of the results of taking the second form is less effective than 
the first one. In Table [Til we give three different kinds of results: without the intrinsic momentum distribution 
(denoted as SI), i.e. S = 1; with the first distribution as Eq. ((28l) (denoted as S2); with the second kind as 
Eq. ([7^ (denoted as S3). Comparing the different results in Table [III we find that at small center mass energy 
y/s = 3.67 GeV the suppression from transverse momentum distribution is more effective: the suppression is 
50% for S2 and 20% for S3. Since the results depend on the explicit form of transverse momentum distribution, 
more experimental results are needed. 

In this calculation, we only present the leading order calculations. The complete next-to-leading order cal- 
culations are much more complicated |34| . For a simple estimate of the size of the next-to- leading order 



contribution, we use the traditional method varying Aqcd and factorization scale t in Eq. ([SO]) . (|49|) and ((53|) : 
Aqcd ~ (0.25 ± 0.05) GeV; changing hard scale t from 0.75t to 1.25t (not changing 1/bi). We find that our 
results are not sensitive to these changes. This implies that the next-to-leading order contribution is probably 
not very large. 



IV. CONCLUSIONS 



In this paper, we have studied the exclusive processes e+e" VP in PQCD approach based on the fcy fac- 
torization. We give three different kinds of results corresponding to different transverse momentum distribution 
functions. With the proper distribution function, our results can be consistent with the experimental results. 
The two different transverse momentum distribution functions S2 and S3 can give about 50% and 20% sup- 
pression respectively at center mass energy -y/s = 3.67GeV. We have included the gluonic contribution for the 
processes involving rj^ ^ meson whose effect is found tiny. We have also included the contribution in which the 
flavor singlet vector meson is produced by an additional photon. This contribution could be neglected at center 
mass energy y/s = 3.67GeV, while these diagrams could induce about 10% difference between e+e^ p^tt" 
and e+e^ p+7r~ at y/s — 10.58 GeV. The s dependence of the cross section has been directly studied which 
indicates that the 1/s'^ scaling is more favored than 1/s^. 
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